arXiv:1502.03334vl [astro-ph.SR] 11 Feb 2015 


Astronomy Astrophysics manuscript no. paper ©ESO 2015 

February 12, 2015 


Colliding Winds in Low-Mass Binary Star Systems: wind 
interactions and implications for habitable planets 

C. P. Johnston^, A. ZhilkidUll, E. Pilat-LohingeiEIH D. BisikalcP, M. GlideP, and S. EggP 

^ University of Vienna, Department of Astrophysics, Tiirkenschanzstrasse 17, 1180 Vienna, Austria 
^ Institute of Astronomy, Russian Academy of Sciences, Moscow, Russia 
® Chelyabinsk State University, Chelyabinsk, Russia 

University of Graz, Institute of Physics, Universitatsplatz 5, 8010 Graz, Austria 
® IMCCE Observatroire de Paris, UPMC, Universite Lille 1, 75014 Paris, France 

February 12, 2015 


ABSTRACT 


Context. In binary star systems, the winds from the two components impact each other, leading to strong shocks and 
regions of enhanced density and temperature. Potentially habitable circumbinary planets must continually be exposed 
to these interactions regions. 

Aims. We study, for the first time, the interactions between winds from low-mass stars in a binary system, to show the 
wind conditions seen by potentially habitable circumbinary planets. 

Methods. We use the advanced 3D numerical hydrodynamic code Nurgush to model the wind interactions of two 
identical winds from two solar mass stars with circular orbits and a binary separation of 0.5 AU. As input into this 
model, we use a ID hydrodynamic simulation of the solar wind, run using the Versatile Advection Code. We derive the 
locations of stable and habitable orbits in this system to explore what wind conditions potentially habitable planets 
will be exposed to during their orbits. 

Results. Our wind interaction simulations result in the formation of two strong shock waves separated by a region of 
enhanced density and temperature. The wind-wind interaction region has a spiral shape due to Coriolis forces generated 
by the orbital motions of the two stars. The stable and habitable zone in this system extends from ~1.4 AU to ~2.4 AU. 
Habitable planets have to pass through strong shock waves several times per orbit and spend a significant amount of 
time embedded in the higher density matter between the shocks. The enhanced density in the wind-wind interaction 
region is likely to lead to a 20% decrease in the size of a planet’s magnetosphere. 

Conclusions. Our results indicate that wind-wind interactions are likely to influence the magnetospheres and upper 
atmospheres of circumbinary planets and could have moderate implications for the development of habitable planetary 
environments. 


1. Introduction 

The formation and evolution of a habitable planetary 
atmosphere is highly sensitive to the surrounding stel¬ 
lar environment. The classical view of habitability is 
based on the ability for a planet to possess liquid water, 
which depends on the amount of visible light from the 
central star incident on the planet ( [Kasting et al.|[T993[ 
Kopparapu et al.]|2013| ). Beyond visible light, stellar out¬ 
put also includes high energy radiation, such as EUV and 
X-rays, high energy particles, and stellar plasma in the 
form of coronal mass ejections and stellar winds. All of 
this output can significantly influence the development 
of habitable planetary environments in ways that are 
currently poorly understood. For example, stellar X-rays 
and EUV can ionise and expand the planetary atmospheres 
and drive hydrodynamic mass-loss ( Watson et al.llTMTl 
ILammer et al.ll?003l ITian et ^12005 ITian et ^120081 
ILammer et al.ll?018|l and winds can compress the plane¬ 
tary magnetospheres and strip away atmospheric particles 
(jKhodachenko et al. 1120121 [Kislyakova et al.]|2014p . 

Since a large fraction of low-mass stars exist in multiple- 
star systems ( |Duquennoy fc Mayor||199H ILadall2006]l . an 
important question is the possibility of planetary habitabil¬ 
ity in binary systems. This is a more complicated question 


due to the gravitational interactions between the planet 
and the two stars and the presence of output from both 
stars. Planets in binary systems can have two types of or¬ 
bital configurations: these are called S-type and P-type. 
S-type orbits are circumstellar orbits where the planet or¬ 
bits one of the stellar components of the binary system. 
P-type orbits are circumbinary orbits where the planet or¬ 
bits both of the stars. Primarily as a result of the Ke¬ 
pler mission, planets in both types of binary systems have 
Doyle et ar]|2011[ IWelsh et al.ll2012l 


been discovered (e.g. 


[Dumusque et air]|2012 


IKostov et al. 11201^ . Although sta¬ 
ble planetary orbits are possible in both types of systems, 
a necessary condition for habitability is that the regions of 
orbital stability co i ncide with the habitable zones. For ex¬ 
ample, Eggl et al. (20131 studied dynamical and radiative 


star-planet interactions in 19 binary systems with differ¬ 
ent stellar masses and orbital parameters and found that 
in two of the systems, dynamically stable orbits are not 
possible within the habitable zones. An additional factor 
that can potentially affect the atmospheres of planets in 
P-type orbits is tidal interactions, which can influence the 
rotational evolution of both stars, and therefore their levels 
of magnetic activity. Since high energy radiation can influ¬ 
ence the evolution of a planetary atmosphere, this could be 
a significant factor in determining the formation of hab- 
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itable environments in a subset of tight binary systems 
(IMason et al.ll20i3ll . A currently unexplored topic is the in¬ 
fluences of winds, and specifically wind-wind interactions, 
on potentially habitable planets in binary systems. This is 
especially interesting for P-type orbits since circumbinary 
planets will pass through wind-wind interaction regions ev¬ 
ery orbit. 

As we discuss in detail in Section colliding winds 
in binary star systems consisting of high-mass stars have 
been researched in great detail. These systems are known 
as colliding-wind binaries (CWBs), and typically contain 
combinations of Wolf-Rayet (W-R) and OB-type stars. 
The fundamental physical mechanisms driving the winds 
of high-mass stars differ significantly from those of low- 
mass stellar winds. Unlike the solar wind, the winds of mas¬ 
sive stars are radiation driven, primarily through line ab¬ 
sorption (jLamers fc Cassinellilll99^ . Such winds are called 
line-driven winds and can have speeds of several thou¬ 
sand km s“^and mass loss rates that exceed 10“® Mq yr“^, 
nine orders of magnitude higher than the mass loss rate of 
the current solar wind. Extreme mass loss of this type can 
have significant effects on the evolution and internal struc¬ 
ture of high-mass stars and can be easily observed (e.g. 
IPuls et aPllQQGI IMokiem et al.llMoT]) . 


The driving mechanisms responsible for the winds of 
low-mass stars are currently poorly understood. Such winds 
are likely to be analogous to the solar wind, which has 
a mass loss rate ~ 10“^^ Mq yr“^ ()Cohenll2011ll and is 
known to be driven by the solar magnetic held. Although 
other acceleration mechanisms, such as wave pressure, are 
probably required, the solar wind is thought to be acceler¬ 
ated primarily by thermal pressure gradients (|Parker II 19581 
ICranmerll2004|l . This requires that the wind is heated as it 
expands into and beyond the solar corona. The main can¬ 
didates for the heating of the wind are energy dissipation 
in the corona by Alfven waves generated by turbulent con¬ 
vective motion in the photosphere and energy dissipation 
by magnetic reconnection of stressed closed magnetic l oops. 
For a review of these mechanisms, see Cranmer (20091. The 
resulting wind streams radially away from the sun in all di¬ 
rections at all times with typical speeds of 400 km s“^and 
800 km s“^. Despite the fact that the wind speed is non¬ 
isotropic and variable in time, the mass flux in the wind 
far from the solar surface is approximately isotropic and 
constant ( |Wang|[r995| |Wang||2010[ |. 


In this paper, we model for the first time the wind-wind 
interactions in a low-mass binary star system. Our chosen 
system consists of two solar mass stars in a tight 0.5 AU 
orbit. We link our results to calculations of the locations of 
stable planetary orbits and habitable zones to show what 
wind conditions a circumbinary planet in the habitable zone 
would see over the course of its orbit. In Section ^ we give 
an introduction to the subject of colliding-winds. In Sec¬ 
tion 1^ we introduce our numerical models for wind forma¬ 
tion and propagation, and for wind-wind interactions. In 
Section]^ we present the results of our wind model. In Sec¬ 
tion we derive the locations of stable and habitable zones 
in our chosen system. Finally, in Section we discuss our 
results and the implications they could have for planetary 
habitability in tight binary systems. 
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Fig. 1. Cartoon showing the main features of wind-wind inter¬ 
actions in a binary system with two identical winds. The wind- 
wind collision (WWC) region takes the form of two shock waves 
(red) separated by a contact discontinuity (blue) and has a spiral 
geometry due to the orbital motions of the two stars. Wind ma¬ 
terial flowing away from one of the stars at supersonic velocities 
impacts the closest shock and is decelerated and thermalised. In 
the post shock gas, the material is then accelerated away from 
the centre of the WWC region. 


2. Colliding Wind Binaries: Observations and 
Theory 

The colliding wind system that we study in this pa¬ 
per is a low-mass counterpart to the well studied collid¬ 
ing wind binaries (CWB^ which are systems contain¬ 
ing two high-mass stars with strong interacting stellar 
winds. Since the vast majority of high-mass stars are in 
binary star systems (e.g. IChini et al.ll2ni2ll . CWBs are 
very com mon. The concept of CWBs was initially intro - 
duced by |Prilutskii & Usov| (| 19761), |Cherepashchuk (|1976|), 
and Cooke et al. (1978|). The early studies of CWBs were 


dedicated mostly to trying to understand the X-ray emis¬ 
sion properties of these systems, without paying atten¬ 
tion to the detailed structure of the wind-wind collision 
(WWC) regions. It was predicted that given the highly su¬ 
personic terminal wind velocities, strong shocks would be 
produced where the winds collide, with post-shock temper¬ 
atures given by 


3 urrip n 

- -V , 

16 kp 


( 1 ) 


where prrip is the average atomic mass and v is 
the pre-shock wind speed. With typical wind speeds 
of u > 1000 km s“^, this will lead to post-shock gas be¬ 
ing heated to temperatures of T > 10 MK, and given the 
high densities in the winds, significant amounts of X-rays 

^ The term ‘colliding stellar winds’, CSWs, is also common in 
the literature. In this paper, we choose not to use either term to 
refer to low-mass binary systems with colliding winds. 
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will be emitted from the WWC reg ions. Observat ional 
snpport for this idea was given by Pollock (19871 and 


Chlebowski & Garmany| (19911 who studied W-R stars and 
O-type stars respectively. They both found higher levels 
of X-ray emission from stars in binary systems than from 
single stars. Since then, many detailed observational and 
theoretical studies of X-ray emission from WWC regions 
in CWBs have been carried out (e.g. IStevens et al.lflMbl 
IGosset et al.ll2009ll . 

Observational detections of wind-wind interactions 
in CWBs in UV, visual, IR, and radio wavelengths 
have also been achieved (e.g. [Koenigsberger fc Au'er||1985[ 
IStevensl|l993l IMadura et al.ll2012ll . Probably the most 
spectacular of these is the direct imaging of the pin- 
wheel nebula coming from WR 104 (ITiithill et al.llT^Ml 
ITuthill et al.ll200'8ll . WR 104 is a W-R star with a close 
OB com panion orbiting wit h a period of approximately 
220 days. Tuthill et al. (20081 presented near-infrared Keck 
images at 11 epochs of an Archimedean spiral emanat¬ 
ing from this system that can clearly be seen to spin 
with the orbits of the two stars. A similar system is 
WR 140, which consists of a W-R star and an OB star on 
a highly eccentric eight year orbit and shows a high level 
of X-ray emission that is clea rly a result wind-wind i nter- 
actions (IWilliams et al. 10.99011 . Dougherty et al. (20051 pre¬ 
sented high resolution interferometric radio observations of 
WR 140 at 23 epochs covering almost a third of the binary 
orbital period. They observed strong radio emission from 
the most dense part of the WWC region between the two 
stars that moves with their orbits. 

Early physical models for the WWC re- 
gion were developed by 
Girard & Willson] (19871, and 
Using simple ram pressure balance 


Huang & Weigert 

'iiUe L hmwn 


(19821, 

(19881. 

Tliey 


arguments, 

were able to predict, for isotropic and steady winds, the 
approximate locations and shapes of the WWC regions, 
assuming that they are infinitely thin and neglecting 
Coriolis forces from the orbital motions of the two stars. 
The latter assumption is only appropriate when the orbital 
velocities of the two stars are much smaller than the wind 
velocities, which can be true for large binary separations 
and fast winds. In this case, the position and shape of the 
WWC region is determined entirely by the ratio 


MiVi 

M2V2 


( 2 ) 


where Mi and Vi are the mass loss rates and pre-shock 
speeds of the two winds respectively. For values of this ratio 
equal to unity, e.g. for two identical winds, the WWC region 
lies in the plane directly between the two stars. As the ratio 
increases, the wind from the first star increasingly domi¬ 
nates the wind from the second star, and the WWC region 
moves towards the second star and becomes increasingly 
bent around it. The WWC region in the line joining the 
two stars occupies the point where the ram pressures of the 
two winds are equaQ Specifically, the ratio of the distance 
between the first star and the WWC region, i?i, and the 
distance between the second star and the WWC region, R 2 , 
is given by 


Ri MiVi 1 

R2 V M2V2 ^ 


( 3 ) 


While instructive, the above model misses most of the 
physics of wind-wind interactions. In reality, the WWC 
region will not be a single infinitely thin layer, but 
will be a thick region of enhanced density and temper¬ 
ature containing a contact discontinuity and separated 
from the quiet winds of the two stars by two strong 
shock waves (ISiscoe fc Heinemannlll97'4ll . This is illus¬ 
trated for two identical winds in Fig. Two impor¬ 
tant physical effects that often need to be taken into ac¬ 
count in CWB models are cooling between the two shock 
waves and Coriolis forces due to the orbital motions of 
the stars. In order to properly take into account all of 
this physics, numerical hydrodynamical models are nec¬ 
essary. The first p roper hydrodynamical models of CWB s 
were published by Luo et al. (19901, Stevens et al. (1992), 
and Myasnikov &: Zhekov J 1993p . These early models ex¬ 
ploited the cylindrical symmetry that exists about the 
line-of-centres between the two stars when the orbital mo¬ 
tions are ignored to run the simulations in 2D (alterna¬ 
tively, 2D simulations run in Cartesian coordinates can 
accurately reproduce a lot of the physics of CWBs, but 
have the disadvantage that the spherically expanding wind 
densities decrease with r~^ instead of with r“^). Models 
in 3D that were able to take into account orbital mo ¬ 
tion s were first presented by|W alder| (19951, Walder (1998), 
and Pittard & Stevens ( 199^| . Since then, with the devel¬ 
opment of more powerful computers and more advanced 
numerical methods, many studies have used numerical hy¬ 
drodynamic and magnetohydrodynamic models to study 
wind-wind interactions in high-mass star systems in increas¬ 
ing detail (e.g. IBisikalo et al. 1120061 ILemaster et al. 1120071 
IPittardir20091 (Lamberts et al.l[2011l IKee et an2014|l . 


Coriolis forces break the cylindrical symmetry about the 
line-of-centres and lead to the WWC region being wrapped 
around the stars and the formation of a pinwheel nebula. 
When the Coriolis force is taken into account, another pa¬ 
rameter becomes important. This is the ratio of the orbital 
velocity to the wind speed. 


I’orb _ 1 G {Ml + M2) 

'i^wind "^wind » d 

where a is the orbital separation of the two stars, Uwind is the 
wind speed, and Mi and M 2 are the masses of the two stars. 
When this value is high, the shape of the wind velocity is 
highly influenced by Coriolis forces (|Lemaster et al.ir2007p . 
This means that orbital motions are more important for 
binary stars on close orbits. For example, two 20 M 0 stars 
on circular orbits with a separation of 0.5 AU have orbital 
velocities of ^ 250 km s“^, which is a factor of a few lower 
than the lower-limits for typical wind speeds. For two solar 
mass stars on circular orbits with a separation of 0.1 AU, 

^ There are in fact three locations where the ram pressures 
could balance: the most important one is in the region where 
the two winds have already approximately reached terminal ve¬ 
locities, and the other two are in the acceleration regions of the 
two winds (see Fig. 2 of IStevens et al.lll992|l . The latter two are 
unstable, with small variations in the wind properties leading 
to one of the winds being completely crushed by the other. In 
this situation, it might be difficult for the crushed wind to ever 
accelerate to high enough speeds to counter the ram pressure 
of the other wind. Therefore, even for isotropic, steady winds, 
and circular orbits, the wind-wind interactions depend on the 
history of the system. 
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the orbital velocities of the stars are ^130 km s“^, which 
is lower than typical solar wind speeds. For the same sys¬ 
tem, with a separation of 5 AU, the orbital velocities are 
^20 km s“^, which is much lower than the expected wind 
speeds. 

Another physical effect that often needs to be consid¬ 
ered is the influence of cooling in the WWC region, which 
is important for determining the separation between the 
two shock waves. For an adiabatic gas, the two shocks are 
relatively far apart. The influence of cooling is to reduce 
the post-shock gas pressure, and therefore the WWC region 
thickness. When cooling is included, the ratio of the cooling 
time in the shock region to the escape time, x = tcool/tesc, 
becomes important. The escape time is approximately the 
time that it takes for the gas to leave the WWC region. 


Stevens et al. (19921 derived a simple scaling law for y. For 
each wind, the importance of cooling can be characterised 
by 


X - 10 


— 18 ^wind^ 


M 


( 5 ) 


where u^ind is the wind speed in km s“^, d is the distance 
between the star and the closest shock in the line-of-centres 
between the two stars in AU, and M is the mass loss rate in 
M0 yr“^. Typically when x ^ 1) cooling due to radiative 
losses in the WWC region is important, and when x 1, 
the wind in the WWC region can be assumed to be adia¬ 
batic. For example, consider a binary system consisting of 
two stars with mass loss rates of 10“® Mq yr“^ and wind 
speeds of 1000 km s“^. For binary separations of 5 AU and 
50 AU, we get values of x of 0.025 and 2.5 respectively, cor¬ 
responding approximately to radiative and adiabatic winds. 
Given the same wind speeds and binary separations, for 
wind mass loss rates of 10“^ Mq yr“^and 10“^"*^ Mq yr“^, 
which are typical for 0-type stars and the Sun, we get val¬ 
ues of X of 2.5 and 2.5 x 10^ respectively. The dependence 
of X on d means that for binary systems on very eccentric 
orbits, the WWC regions can oscillate between radiative 
and adiabatic states, as can be seen in the simulations of 


Pittard (20091. Also, since this parameter can be very dif- 
ferent for the two winds, it is possible that the wind is 
radiative on one side of the contact discontinuity and adi¬ 
abatic on the other side (IStevens et al.lfTM2l) . The reduc¬ 
tion in the thickness of the WWC regions due to radiative 
cooling has important effects on the wind dynamics. For 
very thin WWC regions, non-linear thin-shell instabilities 
can form, leading to highly unstable and variable WWC re¬ 
gions that can take on complex time-dependent geometries 
(IStevens et al.llTMSl [Myasnikov et al.|l998[ IPittardll2009l 
ILamberts et al.|l201ip. 

The development of more advanced physical mod¬ 
els for colliding wind binaries has allowed detailed com- 
parisons between theory and observation. For example. 


Bisikalo et al. (20061 used a colliding wind model to de- 
scribe the outburst of the symbiotic star Z And as the re¬ 
sult of the development of a wind from the white dwarf 
companion that then collides with the wind from the 
donor star. Of particular difficulty is explaining the ob¬ 
served X-ray properties of CWBs, which for a given mas¬ 
sive star system, is a complex function of the orbital 
and wind parameters and the angle at which the system 
is viewed (e.g. iP ittajd fc Stevensll 19971 ISana et ^120041 
IPittard fc ParkiJ2010p . In the context of X-ray emission. 
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an important issue is whether or not the shocks in these 
systems are collisional. If the shocks are non-collisional, the 
postshock electron temperatures in the WWC region will be 
lower than the ion temperatures. If the rate at which elec¬ 
trons are then heated through Coulomb collisions is low, 
the gas might travel a significant distance before temper- 


ature equilibrium is established. Zhekov & Skinner (20001 
used a 2D hydrodynamical wind model to reproduce the ob- 
served X-ray spectra of WR 140 and found better fits when 
separate temperatures were assumed for the electrons and 
the ions in the postshock gas. Similar results were found 
for WR 140 by [Pollock et al^ (|2005| and for WR 147 by 


Zhekov (20071. 


We have given, in this section, a brief and incomplete re¬ 
view of the research into colliding winds in massive binary 
star systems. This is a field that has seen, and continues 
to see, a huge amount of activity. A subject that has not 
yet been addressed in the literature is the colliding winds 
in low-mass star systems. Although there are indirect ob- 
servational constraints on the winds of low-mass stars (e.g. 
Gaidos et al.||2000 Wood et al. 20051, due to our current 


inability to directly observe the winds of such stars, binary 
wind interactions in such systems are less interesting from 
an observational point of view than they are for the much 
more massive winds of high-mass stars. However, they are 
more interesting from the perspective of star-planet inter¬ 
actions and planetary habitability. This is the subject that 
is discussed in the remainder of this paper. 


3. Colliding Wind Model 


In this paper, we simulate the propagation and interaction 
of stellar winds coming from the two components of a tight 
binary system composed of two solar mass stars orbiting 
each other on circular orbits (cbinary = 0) with an orbital 
distance of 0.5 AU. For simplicity, we do not consider the 
effects of binary eccentricity. The orbital period of our sys¬ 
tem is 91 days with each star travelling at approximately 
60 km s“^. Classical tidal interaction theory for stars with 
convective envelopes (|Zahnlll975l IZahnllT977|l suggests that 
the time-scale for the synchronisation of the stellar rotation 
rates with the orbital period for this system is likely to be 
of the order of several hundred Gyrs and the timescale for 
circularisation is likely to be several orders of magnitude 
longer (meaning that our zero eccentricity system repre¬ 
sents only one possible configuration)!^ Our system is there¬ 
fore not tidally locked and the rotation rates for both stars 
will evolve independently and will therefore depend primar¬ 
ily on the age of the system. We assume that each star pos¬ 
sesses an isotropic wind with properties similar to that of 
the slow solar wind, which has a speed of 400 km s“^ at 
1 AU. This is likely to be a reasonable assumption for older 
(> 1 Gyr) systems where the stars have mostly spun-down 
to levels of magnetic activity similar to the present Sun, 
but is unlikely to be reasonable for young rapidly rotating 
stars. 

For our binary wind interaction simulation, we use a 3D 
numerical hydrodynamical code described in Section[3.1| In 


^ We note that an alternative mechanis m for tidal syn chro¬ 
nisation and circularisation presented by Tassoul (19871 and 


Tassoul (19881 could lead to the synchronisation of the stellar 


rotation rates with the orbital period early within the lifetime 
of our system. 
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order for our simulation to be computationally feasible, we 
do not extend our computational domain to the surfaces 
of the two stars, but instead simulate the stars as spheri¬ 
cal masks with radii of 0.1 AU. At the boundaries of these 
masks, we hold the wind properties constant for the en¬ 
tire simulation and we do not include the regions inside 
the masks in the simulation. We assume that the winds 
accelerate to approximately their terminal velocities within 
these masks, and are therefore inserted into the simulations 
at the boundaries of these masks already fully formed. We 
therefore need to estimate the properties of the slow solar 
wind at 0.1 AU, which we do in Section [3?^ 

3.1. Numerical Method 

We run our binary wind simulation using the magnetohy¬ 
drodynamics code Nurgush (jZhilkin fc Bisikaloll2ninp . Nur- 
gush solves the hydrodynamic equations in three dimen¬ 
sions, using a non-inertial coordinate system that rotates 
with the orbital motion of the two stars. We do not con¬ 
sider the influence of magnetic fields on the dynamics of the 
plasma. Therefore, our simulation corresponds to a plasma 
with a very high plasma-/3, where the plasma-/3 is the ratio 
of the plasma thermal pressure to the magnetic pressure. 
Given the high plasma-/3 values in the solar wind far from 
the Sun, this is likely to be a reasonable approximation. 
The set of equations that we solve is 


dp 

at 


+ V • (pv) = 0 , 


( 6 ) 


dv VP 

— + (v V)v =-— -f 2(v X - V0, 


( 7 ) 


ds P 

— -b (v • v)£-f —V • V = 0, 
Ot p 


( 8 ) 


where p is the mass density, v is the velocity, e is the energy, 
P is the pressure, H is the angular velocity of the rotating 
non-inertial frame of reference which we set to the orbital 
angular velocity of the two stars, and (p is the Roche po¬ 
tential. We solve this system of equations using the higher- 
order G udonov-type di ffere ncing scheme, as describe d in 
detail by Zhilkin (20101 and Zhilkin & Bisikalo (20101. 

We solve the equations on a 3U Gartesian grid with 
512, 512, and 50 grid points in the x, y, and z directions re¬ 
spectively. Taking the point {x, y, z) = (0,0,0) as the centre 
of mass of the binary system, our computational domain 
extends from -5 AU to 5 AU in the x and y directions, 
and from -0.25 AU to 0.25 AU in the z direction. The two 
stars are assumed to be located at y = 0 and z = 0, and at 
X = —0.5 AU and a; = 0.5 AU. The stellar winds are gen¬ 
erated by placing spherical masks of radius 0.1 AU around 
the two stars. On the surfaces of these masks, we specify the 
densities, velocities, and temperatures of the winds, which 
we assume to be fully formed. We hold these quantities 
constant in time by resetting them to their initial values 
every timestep. For the initial conditions, we All the rest of 
the computational domain with a low density and temper¬ 
ature gas. For the boundary conditions on each of the six 
outer boundaries, we assume outflow conditions. Therefore, 
a wind is able to flow from the two spherical masks into the 


computational domain, and after a short time, the initial 
conditions are flushed out. 

In reality, by analogy with the solar wind, the expand¬ 
ing wind outside of the WWG region could be heated, and 
inside the WWG region, both heating and cooling could 
be taking place. In simulations of winds from both single 
stars and binary star systems, where detailed knowledge 
of the required heating and radiative cooling is often not 
available, it is common to simulate heating and cooling by 
varying the value of the adiabatic index, 7 . For a monatomic 
ideal gas, assuming a value of 7 of 5/3, without a heating 
function, corresponds to an adiabatic fluid, and assuming 
lower values of 7 leads to heating of the gas as it expands 
and cooling of the gas as it contracts. Lower values of 7 can 
be used in GWB simulations to cool the post-shock plasma 
in WWG region (e.g. IBisikalo et al.l[T997|l . Gooling in the 
WWG region would lead to a decrease in the thickness of 
the region between the two shock waves and the formation 
of thin-shell instabilities, as discussed in Section Using 
Eqn.[^ we estimate that the parameter x that characterises 
the importance of cooling in the WWG region is ^ 5 x 10® 
for both winds in our simulation. We therefore do not ex¬ 
pect radiative cooling to be significant, and set 7 to 5/3. 
This leads to a large simplification of the problem by allow¬ 
ing us to avoid the thin-shell instabilities that are present 
in highly radiative WWG regions. 

3.2. Stellar Boundary Conditions 

In order to generate winds inside the computational do¬ 
main, we need to define the wind parameters (i.e density, 
velocity, and temperature) at the two stellar boundaries. 
It is typical in simulations of the interactions of radiation 
driven winds in high-mass star systems to assume a /3-law 
for the velocity of the wind, such that the velocity is given 
by v{r) = z;oo(l — Ro/r)^, where Rq is some previously de¬ 
fined base of the wind (possibly the stellar surface), v^o is 
the terminal velocity of the wind, and f3 is a, parameter that 
determines how quickly the wind accelerates to its termi¬ 
nal velocity (jPauldrach et al.lll986|l . We choose instead to 
adopt a description of the winds from both stars that is 
based on the solar wind. This is likely to be more appropri¬ 
ate for magnetically driven winds on low-mass stars. 

To calculate the solar wind parameters at 0.1 AU, we 
construct a ID hydrodynamic wind model in a spheri¬ 
cal coordinate system using the Versatile Advection Gode 
(VAG; nWhllToMl HbthllMll- VAG is a freely avail¬ 
able general purpose MHD code that has been used 
extensively in the past for simulations of the solar 
wind dKeppens fc Goedbloed||1999[ [Zieger fc Hansen||2008[ 
I, Jacobs fc Poedtsll2fllip . Bv only simulating the wind in ID 
using an adaptive mesh that is refined close to the so¬ 
lar surface and becomes increasingly coarse far from the 
Sun, we are able to simulate the acceleration and propa¬ 
gation of the wind from the stellar surface to 1 AU. We 
use measured properties of the slow solar wind at 1 AU 
to constrain the free parameters in the model. Our model 
is similar to previous solar wind models in the literature 
(j, Jacobs fc Poedt^l2011l [Kislyakova et al.|2014[ |. 

Since the fundamental driving mechanisms responsi¬ 
ble for the heating and the acceleration of the solar 
wind are not known, solar wind models typically drive 
the wind by assuming a polytropic equation of state, 
P = Kp'^. As discussed in the previous section, setting 
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Fig. 2. Plots showing our one-dimensional hydrodynamical 
model for the slow solar wind that we use to derive the inner 
boundary conditions for our binary wind simulation. From top 
to bottom, the plots show particle number density, velocity, and 
temperature in the wind against distance from the central star. 
The dashed line in the velocity plot shows the sound speed as 
a function of radius. The vertical dotted lines in each plot show 
the location of 0.1 AU, which is where we take the values for the 
inner boundary condition in the binary wind simulation. 



7 < 5/3 means that the wind is implicitly heated as it 
expands. The case of 7 = 1 everywhere corresponds to 
an isothermal Parker wind ()Parkerlll958|l . In order to ac¬ 
curately reproduce the solar wind, 7 must be assumed 
to be a function of distance from the Sun, with the 
wind being almost isothermal close to the Sun, and al¬ 
most adiabatic further out. Close to the solar surface, 
the corona and the wind are well known to be approxi¬ 
mately isothermal, with typical values for 7 being around 
1.1 USteinolfson fc Hundhaus^ll9881IR.oussev et al.llTOOSjl. 


Using Helios 1 data, Totten et al. (19951 measured a value 
of 7 far from the Sun of 1.46. It would therefore make sense 
for us to use this value. However, in our binary wind model, 
we assume a value of 7 of 5/3 in order not to add unphysi¬ 
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cal cooling in the WWC region. We therefore assume in our 
solar wind model that the wind is adiabatic far from the 
Sun and take 7 = 1.1 close to the Sun. We assume that the 
change between these two values happens around 20 R©, 
although this value is to some extent arbitrary given that 
we have no spacecraft measurements of the solar wind tem¬ 
perature structure within ~0.3 AU. In order to ensure a 
smooth transition between the two values of 7 , we vary it 
from the inner value to the outer value linearly between 
15 R© and 25 R©. Since we assume in our binary wind 
simulation that the gas is adiabatic everywhere, and set the 
stellar boundaries at 0.1 AU from the stellar surfaces, there 
is a slight discrepancy in the values of 7 between our solar 
wind simulation and our binary wind simulation. However, 
this is only a small discrepancy that exists in a small part 
of the computational domain, and therefore is unlikely to 
significantly influence the wind dynamics. 

The two free parameters in the model are then the base 
density of the wind, uq, and the base temperature, Tq. In 
our model, the densities of the wind at all radii, and the 
total mass loss rate, vary linearly with the base density. 
However, the wind velocities are only dependent on the 
base temperature. The base temperature is also very im¬ 
portant for determining the densities further from the star 
and therefore the mass loss rates, and the model is very 
sensitive to small changes in this parameter. Since we want 
to reproduce the slow solar wind conditions, we take a value 
of Tq of 2.9 MK, which leads to a wind of 402 km s“^ at 
1 AU. We take a particle density at the base of the wind 
of 2.2 X 10^ cm“^, which gives a proton number density of 
5 cm“^ at 1 AU, consistent with typical slow wind condi¬ 
tions. Our solar wind model based on these parameters is 
shown in Fig. The wind accelerates mostly very close to 
the star, becomes supersonic within 5 R 0 , and by 20 Rq, 
is already travelling at ~350 km s“^. Between 20 R© and 
1 AU, the wind accelerates much slower up to 402 km s“^. 
Due to the way in which we vary 7 with radius, at around 
20 R©, there is an unphysical bump in the temperature and 
sound speed. This bump is small and does not influence 
our results. At 0.1 AU, this model gives a particle density 
of 730 cm“^, a velocity of 340 km s“^, and a temperature 
of 1.47 MK. We use these values as the inner boundary 
conditions for both stars in our binary wind simulation. 

4. Results: Wind Interactions 

In Section we describe our numerical model for the prop¬ 
agation and interaction of stellar winds from two solar mass 
stars with a separation of 0.5 AU. At the beginning of the 
simulation, a wind is made to propagate from the two stel¬ 
lar boundaries, leading to the formation of shocks. These 
shocks propagate outwards, interact with each other in the 
region between the two stars, and propagate to the edges 
of the computational domain where they leave the simula¬ 
tion. After this, the simulation quickly relaxes to a steady 
state, and what remains is the final solution. It takes ap¬ 
proximately 10,000 timesteps for the simulation to come to 
a steady state, which corresponds to 36 days in simulation 
time. We run the simulations to 90 days and see no further 
change in the final result. The time it takes for the simula¬ 
tion to get to a steady state is approximately 1.5 times the 
time for the winds emanating from the stellar boundaries 
to reach the outer boundaries. Since this initial phase has 
no physical significance, we do not discuss it further. 
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Fig. 3. Cuts in the a:j/-plane showing results for the hydrodynamic binary wind simulation for the entire domain {left panel) and 
for the region close to the star {right panel). The background colour and the contour lines show density contours and the arrows 
show the direction and speed of the plasma in the inertial frame of reference. The dashed red line traces the centre of the WWC 
region. 





X (AU) X (AU) X (AU) 

Fig. 4. Plots showing particle number density {left panel), velocity {middle panel), and temperature {right panel) along the 
line-of-centres directly between the two stars (see FiglH. The dashed line in the velocity plot shows the sound speed along the 
line. The vertical dotted lines represent approximately me locations of the two shocks. 


A cut through the computational domain in the or¬ 
bital plane of the two stars is shown in Fig. The back¬ 
ground colour and the contour lines show density, and the 
arrows show the direction of motion of the gas. The form 
of the solution is that of a wind-wind collision (WWC) re¬ 
gion of enhanced density and temperature, surrounded by 
strong shocks, as illustrated in Fig. Outside of the in¬ 
teraction region, the computational domain is filled with 
a quiet wind propagating outwards with an approximately 
constant radial velocity. The winds from the two stars in 
the quiet regions are completely separated from each other 
by the WWC region. The centre of the two shocks on the 
line-of-centres is directly between the two stars, as would be 
expected for winds with equal ram pressures. The dashed 


red line in Fig. traces the line directly between the two 
shocks. Due to the identical speeds of the two winds, no 
clear contact discontinuity can be seen. Across each shock, 
the density increases by about a factor of four, which is 
what we would expect for a strong adiabatic shock. We 
find a density at the centre of the line-of-centres joining 
the two stars of approximately 500 cm“^. As expected, the 
temperature also increases across each shock, leading to 
temperatures as high as 2.5 MK. These temperatures are 
consistent with what is expected in the post-shock gas be¬ 
hind an adiabatic strong shock given the wind velocities 
(using Eqn. [Hand v = 400 km s“^ gives a temperature of 
approximately 2 MK). 
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Fig. 5. Plots showing particle number density {left panel), velocity {middle panel), and temperature {right panel) along the line 
in the centre of the WWC region (i.e., the dashed red line in Fig.l^ as a function of distance along the line. The two lines in each 
plot correspond to the two directions from the centre outwards. The dashed lines in the velocity plot show the sound speed along 
the lines. The increase in wind speed far from the centre of the system is mostly a result of Coriolis forces in the rotating frame 
of reference of our simulation. 


Wind material leaving one of the stars in the direction 
of the other star travels almost 0.2 AU until it reaches the 
nearest shock. Going through the shock, most of the kinetic 
energy in the wind is converted into thermal energy. The 
two winds in the postshock region compress against each 
other and the gas quickly slows to subsonic speeds by the 
time it reaches the point directly between the two stars. 
The wind properties along the line-of-centres between the 
two stars are shown in Fig. Due to the coarseness of the 
grid used in the simulation, not many grid points lie along 
this line. This is likely to be the reason why the profiles 
shown in Fig. ^are not completely symmetric with respect 
to zero, though they are almost symmetric. A region of high 
pressure gas is produced which is then accelerated away 
from this region in the direction parallel to the two shocks. 
In Fig.[^ we show the wind properties along the line directly 
between the two shocks in the interaction region. The gas 
is quickly accelerated by thermal pressure gradients due 
to the high temperatures in the post-shock gas and reaches 
approximately terminal velocity at around 0.5 AU from the 
centre of the computational domain. 

The wind inside the WWC region ends up approxi¬ 
mately equal to the speeds of the gas outside the WWC 
region. Further out, the velocity continues to increase in the 
frame of reference of the two stars due to Coriolis forces. 
The similarity between the speeds of the wind inside and 
outside the WWC region is an inevitable consequence of 
our assumption that the wind is adiabatic everywhere. Gas 
streams out of the stellar boundaries with a certain amount 
of kinetic energy which is then converted into thermal en¬ 
ergy as the gas passes through one of the shocks. The ther¬ 
mal energy is then converted back into kinetic energy as 
the gas is accelerated away from the centre of the WWC 
region. Although most of the thermal energy is reconverted 
into kinetic energy, the temperature in the WWC region re¬ 
mains a factor of a few higher than the temperature in the 
surrounding wind, and the velocity of the wind is always 
slightly lower. 

The influence of the orbital motions of the stars on the 
wind and the WWC region is clearly visible. The Cori¬ 
olis force twists the WWC region into the shape of an 
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Archimedean spiral in the xy-plane. Assuming that the (p 
coordinate is defined such that the line of ^ = 0 points 
along the positive x-axis and the = tt 12 line points along 
the positive y-axis, the Archimedean spiral can be described 
by r « 6.7 — 4.30, where r is the radial distance from the 
centre of the spiral, implying that the spiral makes full 
turns every ~25 AU. In simulations with equal winds from 
both stars and a circular orbit, the shape of the spiral is 
determined by the ratio of the wind velocity to the or¬ 
bital velocity of the two stars. The orbital velocity of the 
two stars is 60 km s“^, giving a ratio for the speed of the 
wind to the orbital speed of ^ 6. The influence of Coriolis 
forces on the dynamics of coll iding winds in binary s ystems 
has been studied in detail by Lemaster et al. (20071. They 
gave simulations with different values of the wind speed 
to orbital speed ratio, and our simulation shows similar 
results to their simulation with a similar ratio (see sim¬ 
ulation C5 in Fig. 9 of their paper). The Coriolis forces 
cause the 3D structure of the shock regions to have a more 
complex shape that depends on the height above the xy- 
plane (|Lemaster et al.ir2007p . Since our computational do¬ 
main only extends up to 0.5 AU in the z-direction, such 
complex shapes cannot be seen. 


In Fig. 1^ we show how the thickness of the shock region 
in the orbital plane of the two stars depends on distance 
from the centre of the spiral. This is calculated by tracing 
the density along circles centred on the centre of the com¬ 
putational domain with successively expanding radii and 
determining along each circle the locations of the shocks. 
We then define the shock thickness as the length along this 
circle between the two shocks. The thickness of the WWC 
region increases with distance from the centre of the com¬ 
putational domain as the winds expand, with distances of 
several AU between the two shocks by the edge of the com¬ 
putational domain. Typically, the WWC region is an eighth 
of the circumference of each circle, such that a cirumbinary 
planet on a circular orbit will spend approximately a quar¬ 
ter of its orbit within the WWC region. 
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Fig. 6. Plot showing the distance between the two shock waves 
as a function of distance along the WWC region. As in Fig. 
the two lines correspond to the two directions from the centre 
of the WWC region that we trace outwards. 
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Fig. 7. Plot showing the maximum orbital eccentricity of a 
planet as a function of semi-major axis in our binary system. 
The colours have the same meanings as in Fig. Due to grav¬ 
itational perturbations, potentially habitable planets around a 
tight binary system periodically generate and lose orbital eccen¬ 
tricity. However, in our binary system, it is unlikely that this 
eccentricity would have a significant influence on the wind con¬ 
ditions seen by a planet in the habitable zone. 


5. Wind Conditions in the Habitable Zone 

In this section, we derive the location of the habitable zone 
in our binary system and combine this with our wind model 
to calculate the wind conditions a habitable circumbinary 
planet would see. In the following discussion, we assume 
that all planets are orbiting in the plane of the binary sys¬ 


tem and in the same direction as the two stars; we do not 
consider retrograde orbits. In our system, the two stars have 
a separation of 0.5 AU, and so circumstellar (S-type) orbits 
would have to be very close to one of the stars to be sta¬ 
ble. Such orbits are not interesting from the point of view 
of planetary habitability, so we therefore only consider cir¬ 
cumbinary (P-type) orbits. 

The determination of the borders of the habitable zones 
in binary star systems has received a large amount of atten- 


20I2||Kane & Hinkel|2013[ 

Eggl et al.|20I31|Kaltenegger & 

Haghighipour|20I3 Haghig 

ripour & Kaltenegger|2013 For- 


gan 2014[ Cuntz 2014 Jaime et al. 20141, especially since 
several planets in such systems have been discovered by 
the Kepler mission. Calculations for the locations of hab¬ 
itable zones are much more complicated when considering 
binary star systems than when considering single star sys¬ 
tems. Firstly, orbital stability must be considered; obviously 
a planet can only be habitable when it has a stable or¬ 
bit, which in binary systems puts limits on the possible 
orbits that a planet can have. When the locations of the 
stable orbits are known, the locations of habitable orbits, 
based on the classical definition of liquid water, must then 
be estimated. The calculations of the classical habitable 
zones are based primarily on considering the flux of stel¬ 
lar radiation onto the planet (i.e. the insolation), but also 
have some dependence on stellar spectral type. From plan¬ 
etary atmosphere models, estimates have been made for 
the lower and upper limits of the insolation that can lead 
to a planet possessing liquid water ( |Kasting et al.| 19931 
[Kopparapu et all2013 l. For a single star, these insolation 
limits can simply be converted into boundaries for the hab¬ 
itable zone, assuming approximately circular planetary or¬ 
bits. In binary systems, the boundaries for the habitable 
zones are much less clear given that the stellar radiation 
flux that a planet receives is necessarily time variable for 
two reasons. Firstly, strong gravitational perturbations in 
binary systems act on a planet orbiting in the habitable 
zone leading to the planet periodically gaining and losing 
orbital eccentricity. This eccentricity causes the planet to 
receive variable amounts of flux from the two stars. Sec¬ 
ondly, even when the planetary orbit is nearly circular, the 
amount of light a planet receives will vary due to changes in 
the spatial configurations of the two stars and the planet. 
The first effect is especially important for planets on cir¬ 
cumstellar (S-type) orbits, and the latter effect is especially 
important for planets on circumbinary (P-type) orbits. 

The first factor that we consider is orbital stability and 
the growth of planetary orbital eccentricity. In a tight bi¬ 
nary system, circumbinary (P-type) planetary orbits with 
small semi-major axes are unlikely to be stable. As a rule 
of thumb, the semi-major axis (measured from the centre 
of mass of the binary system) must be at least twice the 
distance between the two stars, which means that plan¬ 
etary orbits with semi-major axes below 1 AU are not 
possible in our system. Using the method described in 


Pilat-Lohinger & Dvorak (20021, we estimate the locations 
of stable orbits and the magnitude of eccentricity that a 
planet will gain as a function of planetary orbital distance. 
For a given orbital distance, this is done by initialising a sys¬ 
tem of two solar mass stars on circular 0.5 AU orbits and a 
test planet and then integrating the motions of the bodies 
forward in time. At any orbit, we determine orbital stability 
using the Fast Lyapunov Indicator (jFroeschle et al.lFlOOTp . 
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Fig. 8. Figure showing the locations of the stable and habitable zones in our system, as discussed in Section The distances 
are the initial orbital distance of the planet measured from the centre-of-mass of the two stars. The purple zone is the region in 
which stable planetary orbits are not possible. Planets in the red regions are not habitable due to being either too hot or too cold 
for the formation of liquid water. The blue, green, and yellow regions show the locations of the permanently habitable zones, the 
extended habitable zones, and the average habitable zone respectively . The solid and do t ted ve rtical lines show the inner and outer 
limits of the classical habitable zone calculated using the method of Kopparapu et al. (20141 for single stars, assuming a central 
star bolometric luminosity equal to the sums of the bolometric luminosities of the two stars m our binary system. 


and the eccentricity is determined directly from the simula¬ 
tions. We find that the smallest stable planetary orbit has 
a semi-major axis of 1.13 AU. Planetary orbits with semi¬ 
major axes larger than 1.13 AU are likely to be stable. In 
Fig.0 we show the maximum eccentricity that a planet will 
gain as a function of the semi-major axis. The maximum 
eccentricity is slightly dependent on where exactly we start 
the planet in our simulations. For example, a planet started 
along the a;-axis would have a slightly lower maximum ec¬ 
centricity than a planet with the same semi-major axis that 
is started along the y-axis. In Fig. we define the maxi¬ 
mum eccentricity as the average of the values calculated 
from different starting positions of the planet. Clearly, we 
cannot expect that the planet will be permanently on cir¬ 
cular orbits, but will generate some eccentricity. However, 
this eccentricity is not likely to be very large for a planet on 
a stable orbit in our system. For example, the centre of the 
habitable zone is at ~2 AU where the maximum eccentric¬ 
ity is likely to be ^0.03. Therefore, when considering the 
wind conditions that a planet will see, we make the simpli¬ 
fying approximation that all orbits in the habitable zone of 
our system are circular. 

We now calculate the locations of the classical habitable 


zone in our system. Kasting et al. (19931 defined the hab¬ 
itable zone as the region around a star in which planets do 
not exceed the following insolation limits 


i>s>o, 


( 9 ) 


where S is the planetary insolation at the top of the at¬ 
mosphere, and / and O denote upper and lower limits 
on the insolation that can lead to liquid water on the 
planet’s surface. In the case of a G2 V star, these limits 
are I = 1.107 and O = 0.356 in units of the solar constant 
(1360 W m“^). Due to temporal variations in the stellar 
flux incident on the planet, in some orbi ts, the planet will 
leave the habitable zone during its orbit. Eggl et al. (20121 
ex panded the classical concept of habitable zones derived 
by Kasting et al. (19931. Depending on how much of the or¬ 
bit is”wTtImmieseinsoration limits, we distinguish between 
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Permanantely Habitable Zones (PHZ), Extended Habitable 
Zones (EHZ), and Averaged Habitable Zones (AHZ). The 
PHZ is where the planet remains within the insolation lim¬ 
its (i.e. satisfies Eqn. ^ at all times. The EHZ is where 
the planet leaves the habitable zone during its orbit, but 
still spends most of its orbit within the habitable zone. 
For the AHZ, we follow Williams & Pollard] (2002), who 
used global circulation models to show that a planet leav¬ 
ing the habitable zone due to high orbital eccentricity 
migh t not prevent the planet from being habitable. There¬ 
fore, Eggl et al. (20121 define the AHZ as where the time- 
averaged insolation, (j), is in the insolation limits. For more 


details, see Eggl et al. (20121. 


In Fig. we show these different zones for our binary 
system, where the blue, green, and yellow areas represent 
the PHZ, EHZ and AHZ respectively. The red area shows 
the stable but non-habitable zone and the purple area shows 
the dynamically unstable region. The solid and dotted black 
lines define the inner and outer borders of the classical 
habitable zone for this system tha t we calculate using the 
method of Kopparapu et al. (20141 for a single star system 
with a bolometric luminosity equal to the sum of the bolo¬ 
metric luminosities of the two stars. A comparison of the 
two methods shows that the inner border of the classical 
study lies wit hin the non-habit able zone calculated using 
the method ofjEggl et al. 
rdr 


(20121. This indicates the impor¬ 
tance of including orbital dynamics when determining the 


borders of habitable zones in binary star systems. The hab¬ 
itable zone shown in Fig. [^defines the inner and outer parts 
of the AHZ and the EHZ. The differences between these 
borders and the borders of the PHZ can be explained as a 
combination of higher eccentricities for planetary motions 
in these areas and the variations in the incident radiation 
over the planetary orbit. Both effects together cause parts 
of the orbits in these areas to evolve beyond the insolation 
limits. 


It is important when discussing what a planet sees dur¬ 
ing an orbit to be clear how we define a planet’s orbit. When 
discussing our three body system, there are two frames of 
reference that are useful for us to consider: these are the in- 
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Fig. 9. Fignre showing the location of the habitable zone in 
our binary system in comparison to our wind model. The solid 
colours have the same meanings as in Fig. and the contour 
lines show wind density, as in Fig. 
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ertial frame of reference and the frame of reference rotating 
about the centre of mass of the binary system with the orbit 
of the two stars. In the inertial frame of reference, the two 
stars orbit each other with a period of 91 days (0.25 years) 
and a planet on a 2 AU prograde orbit, orbits with a pe¬ 
riod of 730 days (2 years). Therefore, in this frame of refer¬ 
ence, the planet will pass through the WWC region fourteen 
times over the course of its two year orbit. In the rotating 
fame of reference, the two stars are stationary and a planet 
at 2 AU orbits with a period of 104 days (0.29 years) in the 
opposite direction to it’s motion in the inertial reference 
frame. In this frame of reference, the planet passes through 
the WWC region twice over the course of its orbit. In the 
following discussion, we define an ‘orbit’ of a planet as a 
complete 360° orbit in the rotating reference frame of the 
binary wind simulation. The advantage of using the rotat¬ 
ing reference frame is that in this frame, the solution to 
the binary wind problem that we simulate and present in 
Section is not time dependent. Interestingly, the planet 
passes through the WWC region not because of its own 
orbital motion, but because of the much faster orbital mo¬ 
tions of the central stars, which cause the WWC region to 
catch up to and overtake the planet; in fact, the influence of 
the planet’s own orbital motion is to decrease the frequency 
with which the planet passes through the WWC region. 

The main reason that we consider circumbinary orbits 
in this paper is that all planets on such orbits will pass 
through the WWC region. In Fig. we show the average 
number of days per year (i.e. 365 days) that a circumbinary 
planet on a circular orbit spends within the WWC region 
as a function of orbital radius. For a planet with an orbital 
radius of 2 AU, this is typically around 81 days, and this 
becomes longer for planets with larger orbital radii because 



orbital radius (AU) 



orbital radius (AU) 

Fig. 10. Plot showing the average fraction of each or¬ 
bit {upper panel) and the average number of days per year 
{lower panel) that a circumbinary planet on a circular orbit 
would spend in the WWC region as a function of orbital ra¬ 
dius. The dashed vertical lines show the minimum radius that 
the planet can be on a stable orbit. The line colours correspond 
to the colours described in Fig. 


of the thickening of the WWC region with distance from the 
two stars, as shown in Fig. At all orbital radii, at least 
within the computational domain, planets spend slightly 
less than a quarter of each orbit in the WWC region. 

We show in Fig. 0 the wind properties seen by planets 
with orbital radii of 1.5 AU, 2.0 AU, and 2.5 AU. These cor¬ 
respond to orbits near the inner boundary, near the middle, 
and just outside the outer boundary of the habitable zone 
respectively. At all orbital distances, as the planet passes 
into and out of the WWC region, it experiences large dis¬ 
continuous changes in wind density and temperature. The 
wind density increases by approximately a factor of four as 
the planet moves into the WWC region as we expect for 
an adiabatic strong shock. A planet on a 1.5 AU orbit is 
exposed to wind temperatures of approximately 4 x 10^ K 
outside of the WWC region and 2 x 10® K between the two 
shocks. On the other hand, the wind velocity in the WWC 
region is lower than outside of it, but only by a negligi¬ 
ble amount. The clear increase in the wind ram pressure in 
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Fig. 11. Variations in stellar wind parameters seen by circumbinary planets on circular orbits with orbital radii of 1.5 AU (green), 
2.0 AU (blue), and 2.5 AU (red). The line colours represent which zone the planet is likely to be in and have the same meanings 
as in Fig. 


the WWC region is due entirely to the increase in the den¬ 
sity. Outside of the WWC region, the planet is embedded 
in the wind of whichever star is closer. Since the distance 
to the closest star for a planet on a circular orbit changes, 
the planet experiences small changes in the wind conditions 
even when embedded in the quiet wind. 

The elevated wind densities and temperatures in the 
WWC region are likely to have some influence on the plan¬ 
etary upper atmosphere. The increase in the ram pres¬ 
sure will lead to a more compressed magnetosphere. The 
edge of the planetary magnetosphere pointing in the up¬ 
wind direction is approximately where the ram pressure of 
the wind balances the magnetic pressure from the magne¬ 
tosphere. If the only parameter that changes is the wind 
density, p^, then the magnetospheric radius is given by 
Rm oc 1 assuming a dipole field geometry. Therefore, 
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a factor of four increase in the ram pressure will lead to 
a magnetosphere that is 80% of its original size. For more 
extended planetary atmospheres, such as those that could 
be present on younger planets with more massive atmo¬ 
spheres that are exposed to higher levels of stellar EUV 
emission than the current Earth, this can lead to more of 
the atmosphere being exposing to the stellar wind. At the 
same time, higher particle densities in the wind will in¬ 
crease the probability that atoms from an extended plan¬ 
etary atmosphere will interact with the wind, leading to 
the production of ions and energetic neutral atoms (ENAs) 
through charge exchange, and the sweeping away of plane¬ 
tary gas dKislyakova et al.'||2013D . Any circumbinary planet 
will be exposed to these enhanced wind conditions for ap¬ 
proximately a quarter of each orbit, and this may influence 
the mass loss rate and the long-term evolution of the plan- 
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etary atmosphere. However, it is unclear whether this in¬ 
fluence is large enough to have a signiflcant impact on the 
atmosphere of a planet. 

6. Discussion 

In this paper, we have studied the interactions of two iden¬ 
tical winds in a binary star system composed of two solar 
mass stars on circular orbits with a binary separation of 
0.5 AU. To do this, we have used the results of a ID model 
for the slow solar wind as input into a 3D hydrodynami- 
cal wind interaction model including Coriolis forces due to 
the orbital motions of the two stars. We have calculated 
the locations of stable and habitable zones in this system 
to work out what wind conditions circumbinary habitable 
planets see as they orbit the binary system. 

In the slow solar wind model, shown in Fig. the wind 
accelerates to approximately terminal velocity by 0.1 AU, 
meaning that in the our chosen system, the wind interaction 
takes place outside of the acceleration regions of the two 
winds. In binary systems with much smaller separations, it 
is possible that this is not the case, leading to more com¬ 
plex wind interactions. Our wind interaction model results 
in a stable steady-state solution that consists of a wind- 
wind collision region with enhanced density and tempera¬ 
ture, separated from the quiet winds of the individual stars 
by two strong shocks. Given that the orbital velocities of 
the two stars are a signiflcant fraction of the wind velocities, 
orbital motions have a signiflcant effect on the geometries 
of these shocks. In systems with much larger separations, 
and therefore lower orbital velocities, the orbital motions 
would not be so important. 

In this paper, we have chosen for simplicity to study 
one possible low-mass binary star system using a simpli¬ 
fied isotropic stellar wind model. In reality, low-mass bi¬ 
nary systems will have ranges of orbital separations, or¬ 
bital eccentricities, and combinations of stellar masses and 
wind properties. At the same time, the winds themselves 
are likely to be more complicated than we have assumed, 
possessing non-isotropic and time-dependent distributions 
of wind properties. At the same time, transient events, such 
as coronal mass ejections, are likely to additionally compli¬ 
cate the situation of winds in low-mass binary star systems. 
These aspects will be the subject of further study. 

Although we have not yet combined our results with 
models of planetary magnetospheres and atmospheres, it is 
interesting to speculate about how the interactions between 
the two winds in a binary system could influence a circumbi¬ 
nary planet. The most interesting planets are those that 
orbit the binary system in the habitable zone. The habit¬ 
able zone in our system extends from ~1.4 AU to ^2.4 AU 
from the centre of mass of the binary system. Inside of 
~I.l AU, stable planetary orbits are not possible. Due to 
gravitational perturbations, planets on stable orbits in our 
system would not remain on circular orbits, but would in¬ 
stead periodically develop and lose eccentricity, though the 
maximum amount of eccentricity that a planet in the hab¬ 
itable zone would gain is not likely to be signiflcant for the 
wind properties seen by the planet. Combining these re¬ 
sults with the results of the wind-wind interaction model, 
we are able to predict the wind conditions that a planet 
orbiting in the habitable zone would see. This is shown in 
Fig.[TT] A potentially habitable planet would have to pass 
through discontinuous shocks four times per orbit in the 


rotating frame of reference of our binary wind simulations 
(corresponding to 28 times per orbit in the inertial frame 
of reference) and spend almost a quarter of each orbit in 
regions of enhanced wind density and temperature. 

The most obvious influence of the WWC region on a 
planet is due to the increase in the wind ram pressures in 
the dense wind-wind interaction region by about a factor 
of four. This will cause the planetary magnetosphere to be 
compressed to ~80% of its size outside of the WWC region. 
For the current Earth, this is unlikely to lead to increased 
atmospheric loss due to interactions with the wind, given 
that the current Earth’s atmosphere is a thin layer around 
the planetary surface, protected by a much larger magne¬ 
tosphere. However, for a much thicker hydrogen dominated 
atmosphere, such as the one that could have existed on the 
primordial Ea rth (|Hayashi et al.||197^ [Sekiya et al.||1980l 
ILammer et al.ll2014p . this extra magnetosoheric compres¬ 
sion could lead to additional stripping of the outer layers 
of the planetary atmosphere. Therefore, it could be that in 
young binary systems, wind-wind interactions have a sig¬ 
niflcant influence on the initial development of the a plan¬ 
etary atmosphere. This is especially true given that young 
stars show higher levels of magnetic activity and probably 
stronger winds, leading to more magnetospheric compres¬ 
sion and therefore higher levels of atmospheric loss. 

An additional influence on a planetary atmosphere by 
the increase in the wind density could be from an increase 
in the creation of energetic neutral atoms (ENAs) from 
charge exchange betwee n planetary hydro gen atoms and 
high speed wind protons. Chassefiere| (1996) suggested that 
ENAs, which will not be influenced by the planetary mag¬ 
netic field, could enter the atmosphere and cause additional 
heating as they lose their kinetic energy. This heating could 
lead to extra mass loss from the planetary atmosphere. 
An additional influence of the wind-wind interactions on a 
planetary magnetosphere could come from the influence of 
the discontinuous changes in wind properties that a planet 
will see as it passes through the shock waves that confine 
the WWC region. An analogous situation is the influence 
of coronal mass ejections (CMEs) on the Earth. During the 
passage of a CME, the Earth’s magnetosphere experiences 
reconnection events and geomag netic storms that c an be 
measured at the Earth’s surface (Gosling et al. One 

difference between these two situations is that CMEs im¬ 
pact the Earth from the dayside, and therefore upstream, 
direction, whereas the shocks at the edges of the WWC re¬ 
gions in binary systems approach the planet from the direc¬ 
tion in which the planet is moving due to it’s orbital mo¬ 
tion, which is mostly perpendicular to the upstream direc¬ 
tion. Whether or not this difference is signiflcant is unclear 
given that in both cases, the planet sees the discontinuous 
changes in the wind properties from the upstream direction. 

If we were to consider an Earth-mass planet with an at¬ 
mosphere and magnetosphere similar to those of the current 
Earth in the habitable zone around a tight binary system, 
our results suggest that the wind-wind interactions would 
have little influence on the atmosphere. On the other hand, 
the influence on a much more extended hydrogen domi¬ 
nated proto-atmosphere around young magnetically active 
stars could be much larger, due to some of the processes dis¬ 
cussed above. Although the situation is therefore too com¬ 
plicated to make definite predictions about the evolution of 
habitable environments in tight binary systems, we suggest 
that the influence of wind-wind interactions on the evolu- 
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tion of habitable planetary environments is a subject that 
deserves further attention. 
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